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A spin metal-oxide-semiconductor field-effect-transistor (spin MOSFET), which combines a 
Schottky-barrier MOSFET with ferromagnetic source and drain contacts, is a promising device 
for spintronic logic. Previous simulation studies predict that this device should display a very high 
magnetoresistance (MR) ratio (between the cases of parallel and anti-parallel magnetizations) for 
the case of half-metal ferromagnets (HMF). We use the non-equilibrium Green's function (NEGF) 
formalism to describe tunneling and carrier transport in this device and to incorporate spin relax- 
ation at the HMF-semiconductor interfaces. Spin relaxation at interfaces results in non-ideal spin 
injection. Minority spin currents arise and dominate the leakage current for anti-parallel magne- 
tizations. This reduces the MR ratio and sets a practical limit for spin MOSFET performance. 
We found that MR saturates at a lower value for smaller source-to-drain bias. In addition, spin 
relaxation at the detector side is found to be more detrimental to MR than that at the injector side, 
for drain bias less than the energy difference of the minority spin edge and the Fermi level. 

PACS numbers: 72.25.2b, 73.23.2b, 73.40.Sx, 85.75.2d 



I. INTRODUCTION 

In recent years, a vigorous research effort to demon- 
strate spintronic devices [TJ [2] has been pursued. One 
of the motivations has been that spin-based devices are 
identified as one of the most promising alternatives to tra- 
ditional, charge-based logic devices by the International 
Technology Roadmap for Semiconductors [3]. Simula- 
tions have predicted that spintronic logic can scale in its 
size with smaller switching energy and less overall power 
dissipation than electronic logic [4]. 

The concept and operating principles of the first mag- 
netic three-terminal device, i.e. the spin current mod- 
ulator, was proposed by Datta and Das [5] in 1990. It 
comprises of a gate controlling the spin precession in a 
semiconductor channel, a ferromagnetic (FM) source in- 
jecting highly polarized spins, and a FM drain detecting 
the spin polarization. The current depends both on the 
relative directions of magnetization of the source and the 
drain and on the gate bias. Without the gate field, the 
spin current modulator can acts as a giant magnetoresis- 
tance device [S]. The gate exerts an effective magnetic 
field (Rashba field) [8] which causes precession of the 
spins as they move along the channel. With enough chan- 
nel length and strong enough Rashba field, the angle of 
spin precession can be varied from or n, generating 
the device ON and OFF states, respectively, for paral- 
lel magnetizations. Since the device operation involves a 
precise phase of the spins, it is desirable to suppress any 
scattering mechanisms, i.e., the device works best at low 
temperature and in samples with few defects. Physical 
realization of the Datta-Das spin current modulator is 



mainly impeded by the substantial difficulties of obtain- 
ing efficient room-temperature spin injection from the 
FM contacts into the semiconductors like GaAs [9] and 
Si [TUj ■ A recent realization of spin field effect transistor 
(FET) [11] using hot-electron transport through FM thin 
films for all-electrical spin polarized injection and detec- 
tion JT3] with a Si channel is encouraging. There the elec- 
tric field controls the transit time of electrons and thus 
its precession in the magnetic field. However it seems un- 
likely that the Datta-Das spin current modulator or any 
other spin precession devices can provide the ON/ OFF 
current ratio comparable to traditional electronic metal- 
oxide-semiconductor field-effect-transistors (MOSFET). 
Another type of a spin transistor, a spin MOSFET, 
was proposed by Sugahara and Tanaka [13 . It is in 
essence a Schottky barrier (SB) MOSFET, where the 
source and drain are ferromagnetic. A half-metal fer- 
romagnet (HMF) was employed in the original proposal, 
i.e., a material having 100% of electrons with one direc- 
tion of spin at its Fermi level [14] . This property is con- 
ducive to higher spin polarization of injected carriers. Sil- 
icon's mature technology base makes it a preferred choice 
for a channel material. Furthermore, the low spin relax- 
ation rate due to its relatively small spin-orbit effects and 
negligible hyperfine interaction gives propagating elec- 
trons in Si a substantially long spin lifetimes T5J. The 
authors of [13 argued that spin MOSFETs might be used 
for high-density non- volatile memory, whose cell contains 
a single spin transistor, as well as for non- volatile, re- 
configurable logic circuits [16] . The difference between 
the spin MOSFET and a spin current modulator is that 
the former does not rely on the phase of the spin pre- 
cession. Instead, the current is controlled by the height 



of the Schottky barriers, which is different for electrons 
with different spins. This is due to the fact that the 
states of electrons with spin along and opposite to the 
magnetization are split by the value of the exchange in- 
teraction in the ferromagnet. The role of the gate is to 
change the electrical potential and thus the thickness of 
the Schottky barriers. Together the directions of magne- 
tization and the gate bias determine the current through 
the transistor. 

In all types of spin transistors, the key device metric is 
the magnetoresistance (MR) ratio, i.e., the ratio of cur- 
rents for parallel and anti-parallel magnetizations. It is 
a measure of the control of carrier transport by the mag- 
netic state of the device. Theoretically, HMFs will have 
100% spin polarization (all the spins being aligned with 
the magnetization) as first predicted by density func- 
tional theory for NiMnSb [14] and also supported by ex- 
periments such as spin-resolved positron annihilation [17 
and infrared reflectance spectroscopy |18| . Therefore, one 
expects injection of 100% spin polarized carriers from a 
HMF source, which would lead to an extremely large MR 
ratio in the spin MOSFET [13]. The class of Heusler al- 
loys of type X-MnSb are excellent candidates for HFM 
19J. Genuine half-metallic interfaces of NiMnSb with 
III-V semiconductors (e.g. InP and CdS) were predicted 
in the [111] direction with both HMF and semiconduc- 
tor being anion terminated at the interface [50]. Alter- 
natively, one could also consider non-magnetic semicon- 
ducting Heusler alloys (e.g. NiMbSb, NiScSb, NiTiSn 
and CoTiSb) which provide a smaller lattice mismatch 
with Heusler type HMF. Of particular importance is Co- 
TiSb, a semiconductor with indirect bandgap with con- 
duction energy minimum along the six-fold degenerate 
TX symmetry lines, just like Si [21 . 

However, high spin polarization of injected carriers has 
not yet been achieved with HMF. Experiments on spin- 
polarized photocmission |22] and spin-polarized tunnel- 
ing [23 with HMF result in spin polarization values far 
below 100%. This is attributed to the presence of a 'mag- 
netically dead layer', i.e., an area close to the surface of 
the FM which is not ferromagnetic. Its random magne- 
tization interacts with the spins of injected carriers and 
thus decreases their spin polarization. In particular, a 
recent first-principle study of half-metallic NiMnSb/CdS 
interface reveals that the NiMnSb surfaces are not half- 
metallic, even if they are stoichiometric and perfectly or- 
dered |20j . Therefore, in order to estimate a realistic per- 
formance of spin MOSFET, one has to take into account 
spin relaxation processes at these interfaces. 

Many publications have addressed the simulation of 
spin-dependent transport. The drift-diffusion approach 
to spin transport with spin flip processes is reviewed 
in |24j . However, quantum tunneling processes can- 
not be handled by this method. The quantum conduc- 
tance treatment had been applied [25] to simulation of 
a spin-flip transistor. Non-equilibrium Green's function 
(NEGF) treatment of tunneling in ferromagnetic metal- 
oxide multilayers with spin relaxation was described in 



[2B1 HZT] . NEGF has also been applied to spin trans- 
port in carbon nanotubes [28] and molecures [29] . The 
contribution of our work is to treat quantum transport 
through metal and semiconductor structures with spin 
relaxation at the interfaces and under the influence of a 
self-consistent electrostatic field. 

In this paper, we describe a full quantum-mechanical 
model for simulating carrier transport in a spin MOS- 
FET, based on the Keldysh non-equilibrium Green's 
function approach [30, 3T1I32J. In our model, we capture 
the physics of carrier injection and extraction, tunnel- 
ing through Schottky barriers, quantum interference of 
electron wave reflections, and spin relaxation. The influ- 
ence of the 'magnetically dead layer' is incorporated via 
a scattering self-energy of interaction of spin of carriers 
and localized electrons, derived within the self-consistent 
Born approximation [26] [33] . We quantify the effects 
of spin relaxation on the MR ratio of the spin MOS- 
FET. The rest of this paper is organized as follows. In 
section II, we introduce the NEGF formalism with a de- 
tailed mathematical description of the physical quantities 
used in the formalism. In section III, we apply this ap- 
proach to the study of spin MOSFET in the coherent 
regime (without spin relaxation). The bias dependence 
of MR ratio is explained in this section. In section IV, 
we examine the effect of the 'magnetically dead layer' on 
spin transport in spin MOSFET. In section V, we ex- 
plore the dependence of the MR ratio on the strength 
of interaction with a 'magnetically dead layer'. Conclu- 
sions are drawn in Section VI. In Appendix A, we present 
a derivation of the scattering self energy within the self- 
consistent Born approximation. In Appendix B, a simple 
ohmic model analysis of spin MOSFET with HMF and 
FM source/drain contacts is discussed. 



II. MODEL DESCRIPTION 

A schematic drawing of a spin MOSFET is illustrated 
in FigJT] We employed here a double-gate structure with 
a thin film semiconductor for which the gate control of 
electrostatics is optimal. The channel is a semiconductor 
and the source/drain contacts are half-metal ferromag- 
nets (HMF) [2T] with magnetization of Mj and M^, re- 
spectively. The device is large in the y direction. There- 
fore Fig TJdepicts the supercell for our transport problem, 
where the wavefunction solution repeats periodically in 
the y direction. We employ the effective mass approach 
to the description of the bandstructure. The Hamilto- 
nian and the Schrodinger equation reduce to the follow- 
ing Sturm-Liouvillc problem in the longitudinal plane 
(r = (x, z)) (omitting the spin relaxation terms here); 



#* CT (r) 
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V r - [M _1 (r)V r *(r)] 



+ \V(r) + e y + M(r) • a] # ff (r) = etf CT (r) (1) 



where M _1 (r) is the effective mass tensor, with diagonal 
elements m~ l and m~ l [H]. We assume that the trans- 
port mass (m x ), transverse mass (m y ) and quantization 




FIG. 1: (a) A schematic illustration of the three terminal (i.e. 
gate contact Vo> drain contact Vo and source contact Vs) spin- 
MOSFET device. A double-gate structure is employed, (b) A 
zoom-in illustration of the supercell used for the construction of 
device Hamiltonian depicting the various self-energies used in the 
calculation. 



mass (m z ) are spatially uniform within each material. 
M(r) is the magnetization at r, a — {(J x ,a y ,a z ) are the 
Pauli spin matrices and V r = d x i + d z \s.. The eigen- 
value of H is the total energy e, while e y = h 2 k 2 /2m y 
is the transverse energy corresponding to the transverse 
mode solution <f>(y) — exp(ik y y) / \/w . In this work, we 
use the finite difference approach to express our physical 
quantities in matrix representation [30| . Assuming that 
the electrostatic potential is separable, i.e., can be rep- 
resented as V(r) « Vi(x) + V2(z), we can further reduce 
our system Eq. (fTl) to a one-dimensional problem |45j . 

The Keldysh non-equilibrium Green's function 
(NEGF) approach [301 ISU E2] is a powerful technique 
in solving electronic transport in nanostructures having 
open boundary conditions. The infinite problem domain 
(f2) which consists of the semiconductor channel and 
two semi-infinite HMF leads is partitioned into an 
exterior (fl e ) and an interior domain (fij), where only 
the solutions within f2j is to be sought. In this work, 
Cli is composed of the semiconductor channel and a 
small segment of the half-metal ferromagnet at the 
source/drain. In the absence of scattering, the Green's 
function G in f2, is written as, 



G(e») = \(e x )I - H - Vi - E c (e.)] 



-i 



(2) 



T, F + Tir is the contacts self-energy [30 and 
y — €j with £j being the subband energy of the j- 



where S 

ex = e-e^ 

th mode due to the z confinement. For a one-dimensional 
lattice with a nearest neighbor coupling energy of t and 
lattice spacing of a, the contact self-energy is written as, 
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where c; = cos(0i/2), Sj — sin(9i/2), i labels the left 
(source) and right (drain) contacts (i.e., i — L,R), and Qi 
is the magnetization angle with respect to z-axis. We will 
designate the majority spin as 'spin up' and the minority 



spins as 'spin down', kj = [2m x (e — e y — €j — E\)\ 1 / 2 /h 
is the wave- vector in the contact i, and the energy of the 
majority band edge is EJ . A similar identity holds for k\ 
and El . The difference between the energies of the mi- 
nority and the majority spin band edges is the exchange 
splitting in the ferromagnet A s = E^—EJ . The Fermi en- 
ergy (equivalently, the electrochemical potential) in con- 
tact i is designated e l F . The energy bandwidth of occu- 
pied states in HMF is defined as E w = e l F — EJ where Ej 
is the energy of the majority spin band. 

In the ballistic case, states in the device are filled and 
emptied through the contacts. Conventionally, they can 
be defined as the filling and emptying functions (analo- 
gous to the in-scattering and out-scattering functions in 
the case of scattering J3D]), 



5rV)k,e, = [l-/o(e. + e s +e J --ci.)]r < (e.) (4) 

where I\ = i[Sj — St] is the broadening functions of 
the respective contact i = L,R. The electron and hole 
correlation functions are defined as 

G n - p (e x )U y ,^=G(e x )K n '° Ut (e x )U y , e] G(e x y (5) 
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The transverse modes 



can be summed over and we obtain the aggregated elec- 
tron correlation function G n , 

G n (e x ) = G(e x )K n (e x )G(e x y (6) 

with the aggregated filling function defined as 

(7) 



£c n (ez) — /]F [Cx + Cj - e F ) Ti 



where F(e x + ej —e F ) is the Fermi Dirac integral of order 
— \ [33]. The diagonal elements of G n (e x ) are related to 
the charge spectral density at energy e x . Once the total 
charge density is evaluated, the electrostatic potential V\ 
can be obtained using Poisson equation self-consistently. 
To include the effect of spin relaxation, we have to 
modify the Green's function in Eq. (pi. Spin relaxation 
processes arise from the interaction of spins of free car- 
riers with the spins of localized electrons, e.g. in the 
'magnetically dead layer'. The Hciscnbcrg Hamiltonian 
for the spin interaction is 



H, 



Js-S. 



(8) 



where the spin operators for free electrons are s and 
those for localized electrons are S, all in units of h. For 
spin=l/2 these operators are related to the Pauli matri- 
ces s = a/2. The interaction energy is given by J. 

Assuming that the localized electrons are numerous, 
they thus form a reservoir causing an incoherent evolu- 
tion of the free carrier spins. The state of the reservoir is 
described by its density matrix. For the case of spin=l/2 
reservoir, it has the form 



P = 



F u A 
A* F d 



(9) 



where the spin- up and spin-down occupation numbers are 
F u and F d (such that F u +F d = l). 

In this paper, we assume that spin relaxation processes 
are elastic (do not change the energy of free carriers). In 
the self-consistent Born approximation |32j . one can ex- 
press the in- and out-scattering self-energy as a function 
of the electron and hole correlation functions [5S1 [55J , 

S^-(e) = l{e)^ jkl GUe), 

K%(e) = 7(e)*?,- w ^(e), (10) 

where 7(e) is the quantity with the dimension of energy- 



squared proportional to the relaxation rate, which de- 
pends on the number of localized spins and the interac- 
tion energy, $ n ' p are the four-index tensors which pro- 
vides a mapping between electron/hole correlation func- 
tions with the in/out-scattering functions. They are 
products of spin operators and the reservoir density ma- 
trix (JQJ) . Their specific form and derivations are provided 
in Appendix |A"] 



In the case of spin=l/2 reservoir and diagonal density 
matrix (A = 0) the explicit form can be derived 
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The broadening function due to scattering is given by; 

r a (e x )U y ,^ = [sf (e,)|e 9!£ , + sr'Mk,*] (12) 

If we assumed there is no coupling between the different 
e y and tj modes, 

^^)k, ei - 2lr y £ , x _ ex de x % 2 (IS) 

The non-coherent Green's function for a particular trans- 
verse mode (e y ,ej) is given by, 

G(e x )U y ,^ = [(e.)J-ff-Vi-Ec(e.)-E.(e.)|« sl ^]~ 1 (14) 

whereas the electron and hole correlation functions are 
defined as; 

G ' (e x )\e y ,€ j — 

G(e x ) [K n - OU \e x )\ Cv ^+K n ' OUt (^)U v , C] ] G?(e.) t (15) 

The solutions are sought by solving the set of func- 
tions E^e*)^, G"*(e»)| e „ e , and G{e x )\ ey , €j self- 
consistently for each energy e x and modes (e y ,6j). This 
iterative process makes it numerically prohibitive to solv- 
ing realistic transport problems. 

In this work, we shall introduce some simplifications 
to make the numerics more tractable. Firstly, we as- 
sumed that the relaxation rate 7 is energy independent. 
Under the condition where the impurity spin state is un- 
correlated and with equal up and down spin occupation 
probabilty i.e. F u =Fd=\, it can also be shown that the 
function $™=<I>p=<I> [55]. Henceforth, the broadening due 
to scattering is given by; 

r s (€ X )U v , tj = 7$A(ea,) = r s (e x ) 

where A is the local density of state, 

A(e x ) = i [G(e x ) - G(e*) f ] = G"(e.) + G p (e x ) (16) 



We had assumed in Eq. ( 16 1 that the local density of 



states for each transverse modes (e y , tj) is the same and 
is only dependent on the longitudinal energy e x . Simi- 
larly , th e scattering self energy S s (e x ) (computed using 
Eq. ( 13 1) and Green's function G(e x ) are modes indepen- 
dent. The aggregated electron correlation function can 
then be computed self-consistently from, 



G"(e x ) = G(e x ) [K"(e x ) + ££*(£«)] G(c) f 



(17) 



where the aggregated in-scattering self-energy for a spin 
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I impurity is written as, 

tT{e x ) 
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F u G1y(e x ) + s G" T (e :c ) -jGri^a;) 



18) 



The scattering strength of a layer of impurities is de- 
fined by the product 7a, where a is the lattice spacing, 
which physically speaking is the thickness of the inter- 
facial layer in our study. Current is calculated from the 
self-consistent solution of the above equations for any 
terminal i 



F 



I(e x )de x 



where I(e x ) is defined as 



1(C) = \Tt [tT(e x )A(e x ) -F t (e x )G n (e x )] 



(19) 



(20) 



The formalism described in this section allows us to com- 
bine the description of quantum transport in the device 
with the incoherent processes of spin relaxation. 



III. SPIN-MOSFET: COHERENT REGIME 

We consider a spin MOSFET with a double-gated 
structure as depicted in Fig IT] with a 3nm thin film Si 
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FIG. 2: (a) Current- Voltage characteristics for spin-MOSFET 
with the contacts magnetization in parallel configuration plotted 
for Vq = 0.2, 0.4, 0.6V. (b) Current- Voltage characteristics for 
spin-MOSFET with the contacts magnetization in anti-parallel 
configuration plotted for Vq = 0.2, 0.4, 0.6V. 
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FIG. 4: (a) Intensity plot of the majority spin density D u (r) 
(where it is scaled to present color contrast i.e. D u (r) ^) for the 
case where magnetization is in (a) parallel and (b) anti-parallel 
configuration respectively. The terminal bias for both cases are 
Vg = 0.6V and Vd = 0.6V. The energy-resolved current is plotted 
to the left. 
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FIG. 3: MR ratio, MR = (I P - I AP )/I A p, plotted for V G = 
0.2, 0.4, 0.6V. The lines are the least square fitted to data points 
with V D >0.3V for each V G . 



channel. The channel length is Ylnm with gate oxide 
thickness of lnm. Due to the strong body confinement, 
electrons predominantly occupy the doubly degenerate 
valleys along k = (0,0, 1). Its energy dispersion can be 
described within the effective mass approximation with a 
longitudinal mass of 0.91mo (jn z ) and a transverse mass 
of 0.19rao (rn x , m y ). HMFs are employed for the source 
and drain contacts. 

Despite optimistic theoretical projections, engineering 
half metallic interfaces with semiconductors is still in its 
early stages of development [12] . For the purpose of this 
work, we shall assume some reasonable material parame- 
ters for HMF to be used for our theoretical calculations. 
Heusler alloys of type X-MnSb have minority spin bands 
with energy gap ranging from 0.5eV to leV |19j . There- 



fore, we assumed that the HMF's minority spin band 
energy minimum to be OAeV above the metal Fermi 
energy. This value coincides with that of NiMnSb [5D] 
based on a first principle calculation. The metallic na- 
ture of HFM is due to the electronic states occupying a 
large energy bandwidth, E w , which is the energy of the 
conduction band minimum of majority spin band from 
Fermi energy. In this work, we arbitrarily set the metal's 
occupied bandwidth to E'„,=2eV r . We also assumed that 
the majority and minority spin bands in HFM has the 
same transport mass as Si. The choice of different trans- 
port mass in HMF would essentially introduce more re- 
flections at the interfaces, which can also be compensated 
by a larger E w . However, in the general case (e.g. Fe), 
the minority and majority spin bands have to be modeled 
differently 35j. 

The current-voltage characteristics of spin-MOSFETs 
for the case of parallel and anti-parallel magnetiza- 
tion configurations in the coherent regime are plotted 
in FigMa) and (b), respectively. They exhibit simi- 
lar current-voltage characteristics, except that the anti- 
parallel magnetization configurations exhibits a drain off- 
set voltage by an amount of A s -E w . This is due to the 
potential blockade of the majority spin at the drain HMF 
contact in the anti-parallel configuration (see Figfijb)). 
FigH] is an intensity plot of the majority carrier den- 
sity of spin MOSFETs for the case of parallel and anti- 
parallel magnetization configurations. Also shown in 
Fig El is the energy-resolved current (on the left). Os- 
cillations in the energy resolved current for the paral- 
lel case are signatures of tunneling through a barrier, 
e.g. commonly observed in a Fowler Nordheim tunnel- 
ing through Si\Si02\Si sandwiched structure [36]. In 
the anti-parallel case (FigEJb)), the potential barrier at 
the drain HMF contact permits transmission of the ma- 



jority spin only when carrier energy is greater than A s . 
Resonance states due to lateral confinement are observed 
when the carrier energy is less than A s , and these states 
will not contribute any current. These resonance states 
results in the oscillatory behavior in the derivatives of 
the potential profile in the anti-parallel configuration. It 
is also numerically challenging to resolve these states in 
the energy domain due to the relatively fine linewidth in 
these strongly localised resonance levels i.e. a numerical 
challenge faced similarly in modeling of resonant tunnel- 
ing diode. A fine gridding is employed to resolve these 
resonance states so that the charge density can be more 
reliably computed. 

The magnetocurrent ratio (MR) is defined to be 
M R=(I p— I ap) 1 1 ap ■ It serves to quantify the magneto- 
resistance difference between the parallel and anti- 
parallel magnetization configuration states of the device. 
Since the device in the anti-parallel state has an apparent 
drain offset voltage of A s -E w , Iap will generally be less 
than I P . Fig(3]plots the MR of the spin-MOSFET under 
different biasing conditions. The MR follow an approx- 
imately exponential relationship with Vp> at a given Vq 
bias. We can also see that the MR begins to diminished 
(A/_R<10%) when the Vd in the anti-parallel configu- 
ration approaches the current saturation condition i.e. 
Vd = V sat - Consequently, there is a rightward shifts of 
the MR versus Vd curve when the gate bias Vg increases, 
since V sa t increases with increasing Vg- 
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FIG. 5: (a) Spin polarization P in the channel as a function of the 
interface scattering strength characterized by "fa, for spin injection 
through a Schottky barrier in a HMF/Si structure, (b) Current- 
voltage characteristics for spin MOSFET with the contacts magne- 
tization in parallel configuration plotted for Vg=0.6V . All spin cur- 
rents are measured at the drain, (c) Current-voltage characteristics 
for spin-MOSFET with the contacts magnetization in anti-parallel 
configuration plotted for Vq=0.QV . For both cases, the spin relax- 
ation strength at the HMF/scmiconductor interfaces on both the 
detector and injector sides are characterized by -ya = 0.5eV 2 nm. 
The minority spin current is plotted as dashed curve. 



IV. SPIN-MOSFET: INCOHERENT REGIME 

Spin exchange scattering processes between the tun- 
neling electron and the localized spin impurities at the 
HMF/Si interfaces are responsible for the incoherent na- 
ture of the electron transport. These localized spin impu- 
rities lead to decoherence of the electronic spins states. In 
our model, we assume that there are external reservoirs 
constantly maintaining these localized spin impurities in 
a state of equilibrium with a random polarization of spin. 
For a spin=l/2 case, the completely unpolarized states 
corresponds to the density matrix with ^=^=0.5 and 
A=0 [3p], 0B]. 

As discussed in Sec II, coupling between the number 
of available electrons/holes ([G"]/[G P ]) for a particular 
state, and the in/out-flow ([E^]/[S° ut ]) to/from that 
state is related through a four-index scattering tensors 
$. The scattering strength, 7a, is assumed to be energy 
independent. It is proportional to nj, the impurity con- 
centration and to (J 2 ), the averaged interaction energy 
of the impurity layer (see Appendix lAj) . These parame- 
ters can be mapped to available experimental data, and 
they are the only parameters used to characterize the 
spin relaxation of the impurity layer. FiglMa) shows the 
spin polarization P in the channel due to spin injection 
from a HMF contact into a long-channel semiconductor. 
P is defined as; 



P = 






(21) 



Our model is capable of capturing the non-ideal spin 
injection of a HMF due to the presence of a magnetically 
dead layer, by tuning the value of 7a (see Fig[5k). The 
spin polarization in the channel decreases with the 
increase of spin relaxation strength in an approximately 
linear fashion. 

FigjEb and [5fc plots the drain current vs. drain 
voltages at Vg—0-6V for both the parallel and anti- 
parallel configuration, in the presence of spin relaxation 
at the HMF/semiconductor interfaces on both the detec- 
tor and injector sides, characterized by spin relaxation 
strength of 7a = 0.5eV 2 nm. In contrast to the current- 
voltage characteristics in the coherent case (see Figp|, 
the minority spin can now contribute to spin current, 
which significantly modifies the current-voltage charac- 
teristics in the anti-parallel configuration. This 'leak- 
age current' is facilitated through spin relaxation at the 
HMF/semiconductor interfaces. Spin relaxation renders 
the blocking (due to the potential barrier at the drain 
HMF) of the majority spin transmission ineffective. Ma- 
jority carriers injected from the source can now undergo 
spin relaxation at the HMF/semiconductor interface at 
the drain side and become a minority spin. Therefore, 
substantial amount of minority spin current is registered 
in the anti-parallel configuration even when the drain bias 
is less than A, -E,„. 




FIG. 6: Intensity plot of the (a) majority spin density n u (f) and (b) minority spin density n^lf) (where it is scaled to present color 
contrast i.e. n u (r) ' 8 and 10 X ri^r) ' 8 respectively) for the case where magnetization is in parallel. Similar plots for the majority and 
minority spin density for the anti-parallel configuration in (c) and (d) respectively. The terminal bias for both cases are Vq=Vd=0.6V . 
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FIG. 7: Spin current spectral density vs. carrier energy, for the 
parallel directions of magnetization of source and drain. Total cur- 
rent - black dashed line, majority spin - red solid line, minority 
spin - blue shaded area. Top plot - at the source, middle- plot - in 
the center of the channel, bottom plot - at the drain. 



Fig)6[a)-(d) plots the majority and minority spin den- 
sities in the parallel and anti-parallel configuration at 
Vg=Vd=0-QV. Their respective energy resolved current 
is plotted in Fig(7](8] In this set of calculations, the mi- 
nority spin density is approximately an order of magni- 
tude smaller than the majority spin density. In the in- 
tensity plot for the minority spin density we observe sig- 
natures of spin flip process at the spin injector interface 
as a layer of high intensity minority spins. Majority spin 
was scattered into the minority spin states in the injec- 
tor as evanescent states with finite probability to tunnel 
through the Schottky barrier into the channel. There- 
fore, HMF will not have perfect spin injection efficiency 
when spin relaxation mechanisms are incorporated into 



FIG. 8: Same as in Fig. [7] for the anti-parallel directions of mag- 
netization of source and drain. 



the model. 

Fig. [7] and [8] plots the corresponding energy-resolved 
spin current density in each of the respective device re- 
gions i.e. (a) source, (b) channel and (c) drain, when 
the device is in the parallel and anti-parallel magnetiza- 
tion state respectively. Although the spin current is not 
conserved across the device, the charge current at each 
energy is conserved i.e. a property ensured by virtue 
of the self-consistent Born framework. From the energy 
resolved spin current characteristics, we can make the 
observation that more minority spin current is produced 
in the anti-parallel configuration than its parallel coun- 
terpart. The energy resolved minority spin current in 
the anti-parallel configuration is characterized by multi- 
ple resonance peaks. This is in part due to the laterally 
confined majority spins in the channel where the minority 
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is evident from FigMb) that spin relaxation at the 
HMF /semiconductor interfaces at the detector side is 
more detrimental to the MR ratio than that due to the 
spin relaxation at the injector side, for drain bias less 
than A s -E w . For each of these cases, Fig 10 shows the 
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FIG. 9: (a) MR ratio (Ip — Iap)/Iap versus drain biases for Vg = 
0.6V. The spin relaxation strength at the HMF/scmiconductor 
interfaces on both the detector and injector sides are character- 
ized by 7a=0, 0.1, 0.5, 1.0eV 2 nm. (b) MR ratio versus drain bi- 
ases for Vg = 0.6V for spin relaxation at both detector and in- 
jector HMF/semiconductor interfaces (dashed line), at detector 
HMF/semiconductor interface only (solid line with square symbol) 
and at injector HMF/semiconductor interface only (solid line with 
circle symbol). 



spin current was derived from it through spin relaxation 
processes. 



V. INFLUENCE OF SPIN RELAXATION ON 
MR RATIO 

In this section, we discuss MR ratio in the presence of 
spin relaxation. FigMa) shows the MR ratio versus the 
drain voltage bias at Vg=0-6V in the presence of spin re- 
laxation at the HMF/semiconductor interfaces on both 
the detector and injector sides, characterized by spin re- 
laxation strength of 7a = — leV 2 nm. As evident in 
FigMa), the increasing spin relaxation strength at the 
HMF/semiconductor interfaces results in decreasing MR 
ratio over the whole drain voltages sweep, especially for 
drain biases less than A s -E w . Spin relaxation process en- 
hance the probability of conduction through the minority 
spin channel which consequently smears the distinction 
between transport in the parallel and anti-parallel config- 
uration, resulting in diminishing MR ratio. Our analysis 
in this work focuses only on HMF contacts. The impact 
of interfacial spin dephasing on normal ferromagnetic 
contacts could have intrinsically different dependence on 
the spin dephasing processes. In the ohmic regime, one 
could model the interfacial spin dephasing process with a 
resistor element (r~) between the majority and minority 
spin channels at each contact interfaces. A simple elec- 
trical analysis of such a setup for spin-MOSFETs with 
either HFM or normal ferromagnetic contacts shows that 
they have different dependence on r 7 (see Appendix pi). 

Lastly, we studied the impact of spin relaxation at each 
of HMF/semiconductor interfaces. FigMb) plots the 
MR ratio for spin relaxation at the HMF/semiconductor 
interfaces at either the injector or detector side. It 



minority spin current in the parallel and anti-parallel con- 
figurations at different drain biases. As discussed pre- 
viously, the fractional contribution of minority spin cur- 
rent in the parallel case is relatively small and arises only 
when the drain bias is larger than A s -E w . More minor- 
ity spin current is produced in the anti-parallel case in 
this regime. From Fig]10[b), it is evident that more mi- 
nority spin current is produced when spin relaxation at 
the HMF/semiconductor interfaces is at the detector side 
compared to when it is at the injector side. The conse- 
quence is that spin relaxation at the HMF/semiconductor 
interfaces at the detector side is more detrimental to the 
MR ratio, at least for drain bias less than A s -E w . A 
plausible explanation to why spin relaxation at the de- 
tector side produces more minority spin current in anti- 
parallel case could be understood from FigJ6^d). When 
Vd <A s -E w , majority carrier undergoing spin relaxation 
at the detector interface will either be admitted into the 
drain as a minority spin or be reflected back into the 
channel as majority spin. Waves that backscatters at the 
detector interface will then be reflected back and under- 
goes another scattering event. This multiple scattering 
processes enhances the spin relaxation rate. This study 
reveals that proper HMF/semiconductor interface treat- 
ment at the detector side is more pertinent to achieving 
high MR ratio in spin MOSFET for operation regime 
where biasing conditions is less than A s -E w . 

We note that a spin MOSFET uses magnetoresistance 
caused by the change in spin transport through the de- 
vice. However in experimental relaizations, a large role 
is played by fringe magnetic fields near the ferromagnets 
which might induce local Hall effect in the semiconductor 
channel at vicinity of the interface. This results in a spu- 
rious contribution to the magnetoresistance, as it was 
explained in [37j [38j [39]. Measurements based on the 
precession of spins in the channel and the correspond- 
ing Hanle effect provide a rigorous method to detect the 
spin transport [TTJ [T5J [3U]. Practically, one should be 
able to minimize this spurious Hall effect via a proper 
geometrical setup of the spin injection device [41] . 



VI. SUMMARY 

In this paper, we have simulated the operation of a 
spin MOSFET using the non-equilibrium Green's func- 
tion approach. A large spin-splitting energy A s for the 
HMF contacts is beneficial for achieving a large MR ra- 
tio. This is because the spin-splitting energy presents a 
potential barrier at the drain side in the anti-parallel con- 
figuration and blocks the transmission of majority spins, 
resulting in a larger drain 'threshold voltage' of A s -E w 
for the device in the anti-parallel configuration. We high- 
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FIG. 10: Minority spin current versus drain bias for Vq=0. 6V at 
(a) parallel and (b) anti-parallel configuration. The spin relaxation 
strength at the HMF/semiconductor interfaces on either/both the 
detector and injector sides are characterized by -ya = 0.5eV 2 nm. 
We consider the situation where spin relaxation are at both detector 
and injector HMF/semiconductor interfaces (solid line), at detector 
HMF/semiconductor interface only (dotted line) and at injector 
HMF/semiconductor interface only (dash-dotted line). 



lighted and explained the bias dependence of MR ratio of 
spin MOSFET. Next, we demonstrated the incorporation 
of spin relaxation at the HMF/semiconductor interfaces 
via scattering self-energies. Once spin relaxation is in- 
cluded in the model, minority spin current arises and it 
dominates the 'leakage current' in the anti-parallel con- 
figuration even when the drain bias is less than A s -E w . 
This substantially reduces the spin MOSFET MR ratio. 
Lastly, we studied the impact of the 'magnetically dead 
layer' at the injector and at the detector sides in isola- 
tion. It was found that spin relaxation caused by the 
'magnetically dead layer' at detector side is more detri- 
mental to the MR ratio. This can be attributed to the 
fact that spin relaxation due to the 'magnetically dead 
layer' at the detector side is more efficient in flipping the 
majority spin to minority spin current. 
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APPENDIX A: DERIVATION OF SCATTERING 
SPIN TENSORS 

In this appendix we provide a simplified derivation and 
the explicit form in a parti cul ar case of spin=l/2 of the 
scattering tensors in Eqs. (111. It follows earlier papers 
[251 133] i but is presented here for completeness sake. In 



general, the scattering tensor is determined by the Hamil- 
tonian of interaction with the reservoir Hi as follows |33j 



Ni 



7*? iH = J2 Tr ^ H ^ H 



a j j 



(Al) 



where p the density matrix of the reservoir spin, like in 
Eq. (pi, and Nj is the number of reservoir modes. In our 
case, the Hamiltonian & corresponds to the spin-spin 
interaction. The spin-dependent four-index tensors are 
thus 



^ijkl 



Tr 



pS a S ] 



S lj S ik> 



(A2) 



where indices a and (3 run over the projections of the op- 
erators on Cartesian axes, and summation over repeating 
index is implied. 

And the prefactor contains the interaction energy J 
averaged over the spectral range of reservoir modes reso- 
nant with the transition (this averaging is designated by 
angled brackets), which according to [26 is, 



E^ 2 



oc ni 



(J 2 ), 



(A3) 



where ni is the area density of impurities in the layer and 
averaging is done with the account of geometry of the 2D 
layer. 

For the particular case of the spins of localized elec- 
trons (reservoir) equal to 1/2, more explicit expressions 
can be obtained. It is convenient to use the raising and 
lowering Pauli matrices 



= U° X + ^ 



The interaction Hamiltonian can be re-written as 



Hi = J 



V S~ 



+ 



-s 4 



a z S z 



(A4) 
(A5) 

(A6) 



It is easy to evaluate the averages of the spin operators 
in the state of the reservoir pJF. The scattering tensor 
becomes 



®ijki = (F u + F d )sijSi h + Asjs^ fe - Astjsf k 



+A* 



s lj s ik 



^ S lj S ik + FdSijS ik + F u SijS ik , 



(A7) 



and similarly for <1> P . 



APPENDIX B: ELECTRICAL ANALYSIS OF 
HMF/FM CONTACTS IN OHMIC REGIME 

The magnetocurrent ratio is an important device 
metric which provides a measure of the degree of 
distinguishability between the parallel and anti-parallel 
states of the spin transistor through their measured 
current. The Julliere's description of the magnetocur- 
rent ratio in terms of the available effective tunneling 
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FIG. 11: (a) and (b) plot the equivalent circuit for the parallel and 
anti-parallel configuration of a spin MOSFET in the Ohmic regime 
where r a /a is the resistance of the majority/minority spins in the HFM 
contacts, r c the channel resistance and r 7 is the resistance related to 
spin dephasing. The transistor switch in the equivalent circuit model 
has a threshold voltage of A s -E w . 



density of states proves to be inadequate [12] , especially 
when there is interfacial spin flip scattering processes. 
A simple analysis using the electrical circuit equiva- 
lent in Fig 11 'a) and (b) can offer more reliable insights. 



Fig 11 a) and (b) plot the equivalent circuit for the 



spin dephasing. The transistor in the circuit is to model 
the effect of the minority /majority spin blocking due 
to the spin exchange energy at the drain HMF in the 
parallel/ anti-parallel configuration. A simple Ohmic 
model can adequately explains the late turn on of the 
minority current in parallel configuration at Vd=^ s -E w 



as depicted in Fig(5j Recall that in the parallel state, no 
minority spin current can be detected at the HMF drain 
when Ve,<A s -E w , as the source injected states (with 
energy lesser than the source Fermi energy) admitted 
into the drain will cither decay cvancsccntly or revealed 
itself as resonant levels in the channel (see FigjHFd)). In 
addition, it also explains why the minority spin current 
in the anti-parallel configuration can leak effortlessly 
into the drain once the spin dephasing process opens up 
the minority spin channel through g 7 . 

Consider the regime where the drain bias is less than 
A s -E w , the respective resistance for the parallel and anti- 
parallel configuration are written as; 



r P - 2r a + 



1 r c (r c + 2r 7 ) 

2 r c + r~, 



TAP 



2r a + - (r c + r 7 ) 



(Bl) 



which both converges to 2r a + r c /2 when the scattering 
is large i.e. r 7 <Cr c , yielding a zero magnetocurrent ra- 
tio. Maximum magnetocurrent ratio is achieved when 
the scattering is minimal i.e. r 1 ^>r c , with a large mag- 
netocurrent ratio of the order r 1 /r c . We can also derive 
similar expressions for normal ferromagnetic contacts; 



rp 



parallel and anti-parallel configuration in the Ohmic 

regime where r a /p {g a /p) is the resistance (conductance) 

of the majority/minority spins in the ferromagnetic tap 

contacts, r c (g c ) the channel resistance (conductance) 

and r 7 (g 7 ) is the resistance (conductance) related to where i 



2r a rp 



+ 



{2r 1 r a + r T c){2r 1 rp + r T r c ) 
r T (2r~ f r / 3 + r T r c ) + r T {2r 1 r a + r T r c ) 



2r a rp r-y(r a + rp) + r T r c 



+ 



2/v 



(B2) 



-rp- 



By making the assumption that 



the r c <^r^, we arrive at the result that the magnetocur- 
rent ratio for the HMF and FM case varies with the spin 
dephasing conductance g 1 according to g^ 1 and g~ 2 
under these limiting conditions respectively. This sim- 
ple analysis illustrates that the impact of interfacial spin 
relaxation might have a different impact on the FM case. 
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